Search for: All records

Abstract Our understanding of the assembly timeline of the Milky Way has been transforming along with the dramatic increase in astrometric and spectroscopic data available over the past several years. Many substructures in chemo-dynamical space have been discovered and identified as the remnants of various galactic mergers. To investigate the timeline of these mergers, we select main-sequence turnoff and subgiant stars (MSTOs) from the H3 survey, finding members in seven metal-poor components of the halo: Gaia-Sausage/Enceladus (GSE), the Helmi Streams, Thamnos, Sequoia, Wukong/LMS-1, Arjuna, and I’itoi. We also select out a metal-poor in situ population to facilitate comparison to the evolution of the Milky Way itself at these early epochs. We fit individual isochrone ages to the MSTOs in each of these substructures and use the resulting age distributions to infer simple star formation histories (SFHs). For GSE, we resolve an extended SFH that truncates ≈10 Gyr ago, as well as a clear age–metallicity relation. From this age distribution and measured SFH we infer that GSE merged with the Milky Way at a time 9.5–10.2 Gyr ago, in agreement with previous estimates. We infer that the other mergers occurred at various times ranging from 9 to 13 Gyr ago, and that the metal-poor in situ Galaxy built up within only a few billion years. These results reinforce the emerging picture that both the disk and halo of the Milky Way experienced a rapid assembly. 

Abstract We presentAugustus, a catalog of distance, extinction, and stellar parameter estimates for 170 million stars from 14 mag <r< 20 mag and with ∣b∣ > 10° drawing on a combination of optical to near-infrared photometry from Pan-STARRS, 2MASS, UKIDSS, and unWISE along with parallax measurements from Gaia DR2 and 3D dust extinction maps. After applying quality cuts, we find 125 million objects have “high-quality” posteriors with statistical distance uncertainties of ≲10% for objects with well-constrained stellar types. This is a substantial improvement over the distance estimates derived from Gaia parallaxes alone and in line with the recent results from Anders et al. We find the fits are able to reproduce the dereddened Gaia color–magnitude diagram accurately, which serves as a useful consistency check of our results. We show that we are able to detect large, kinematically coherent substructures in our data clearly relative to the input priors, including the Monoceros Ring and the Sagittarius Stream, attesting to the quality of the catalog. Our results are publicly available at doi:10.7910/DVN/WYMSXV. An accompanying interactive visualization can be found athttp://allsky.s3-website.us-east-2.amazonaws.com. 

ABSTRACT Post-common envelope binaries (PCEBs) containing a white dwarf (WD) and a main-sequence (MS) star can constrain the physics of common envelope evolution and calibrate binary evolution models. Most PCEBs studied to date have short orbital periods (Porb ≲ 1 d), implying relatively inefficient harnessing of binaries’ orbital energy for envelope expulsion. Here, we present follow-up observations of five binaries from 3rd data release of Gaia mission containing solar-type MS stars and probable ultramassive WDs ($$M\gtrsim 1.2\ {\rm M}_{\odot}$$) with significantly wider orbits than previously known PCEBs, Porb = 18–49 d. The WD masses are much higher than expected for systems formed via stable mass transfer at these periods, and their near-circular orbits suggest partial tidal circularization when the WD progenitors were giants. These properties strongly suggest that the binaries are PCEBs. Forming PCEBs at such wide separations requires highly efficient envelope ejection, and we find that the observed periods can only be explained if a significant fraction of the energy released when the envelope recombines goes into ejecting it. Our one-dimensional stellar models including recombination energy confirm prior predictions that a wide range of PCEB orbital periods, extending up to months or years, can potentially result from Roche lobe overflow of a luminous asymptotic giant branch (AGB) star. This evolutionary scenario may also explain the formation of several wide WD + MS binaries discovered via self-lensing, as well as a significant fraction of post-AGB binaries and barium stars. 

ABSTRACT We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding time-scales of τin = 1.01 ± 0.13 Gyr for GSE and $$\tau _\text{in} = 3.08^{+3.19}_{-1.16}$$ Gyr for Wukong/LMS-1. GSE formed stars for $$\tau _\text{tot} = 5.40^{+0.32}_{-0.31}$$ Gyr, sustaining star formation for ∼1.5–2 Gyr after its first infall into the Milky Way ∼10 Gyr ago. Our fit suggests that star formation lasted for $$\tau _\text{tot} = 3.36^{+0.55}_{-0.47}$$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass M⋆ predicted by simulations and semi-analytic models of galaxy formation. Our inferred values of the outflow mass-loading factor reasonably match $$\eta \propto M_\star ^{-1/3}$$ as predicted by galactic wind models. Our fitting method is based only on Poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes (20 ≤ N ≤ 2000) and measurement uncertainties (0.01 ≤ σ[α/Fe], σ[Fe/H] ≤ 0.5; $$0.02 \le \sigma _{\log _{10}(\text{age})} \le 1$$). Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space. 

Abstract The Magellanic Stream (MS)—an enormous ribbon of gas spanning 140° of the southern sky trailing the Magellanic Clouds—has been exquisitely mapped in the five decades since its discovery. However, despite concerted efforts, no stellar counterpart to the MS has been conclusively identified. This stellar stream would reveal the distance and 6D kinematics of the MS, constraining its formation and the past orbital history of the Clouds. We have been conducting a spectroscopic survey of the most distant and luminous red giant stars in the Galactic outskirts. From this data set, we have discovered a prominent population of 13 stars matching the extreme angular momentum of the Clouds, spanning up to 100° along the MS at distances of 60–120 kpc. Furthermore, these kinematically selected stars lie along an [α/Fe]-deficient track in chemical space from −2.5 < [Fe/H] <− 0.5, consistent with their formation in the Clouds themselves. We identify these stars as high-confidence members of the Magellanic Stellar Stream. Half of these stars are metal-rich and closely follow the gaseous MS, whereas the other half are more scattered and metal-poor. We argue that the metal-rich stream is the recently formed tidal counterpart to the MS, and we speculate that the metal-poor population was thrown out of the SMC outskirts during an earlier interaction between the Clouds. The Magellanic Stellar Stream provides a strong set of constraints—distances, 6D kinematics, and birth locations—that will guide future simulations toward unveiling the detailed history of the Clouds. 

Abstract Phosphorus (P) is a critical element for life on Earth, yet the cosmic production sites of P are relatively uncertain. To understand how P has evolved in the solar neighborhood, we measured abundances for 163 FGK stars over a range of –1.09 < [Fe/H] < 0.47 using observations from the Habitable-zone Planet Finder instrument on the Hobby–Eberly Telescope. Atmospheric parameters were calculated by fitting a combination of astrometry, photometry, and Fe I line equivalent widths. Phosphorus abundances were measured by matching synthetic spectra to a P I feature at 10529.52 Å. Our [P/Fe] ratios show that chemical evolution models generally underpredict P over the observed metallicity range. Additionally, we find that the [P/Fe] differs by ∼0.1 dex between thin disk and thick disk stars that were identified with kinematics. The P abundances were compared withα-elements, iron-peak, odd-Z, and s-process elements, and we found that the evolution of P in the disk most strongly resembles that of theα-elements. We also find that molar P/C and N/C ratios for our sample match the scatter seen from other abundance studies. Finally, we measure a [P/Fe] = 0.09 ± 0.1 ratio in one low-αhalo star and probable Gaia–Sausage–Enceladus member, an abundance ratio ∼0.3–0.5 dex lower than the other Milky Way disk and halo stars at similar metallicities. Overall, we find that P is likely most significantly produced by massive stars in core-collapse supernovae, based on the largest P abundance survey to date. 

Abstract The majority of the Milky Way’s stellar halo consists of debris from our galaxy’s last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from the GSE have been kinematically and chemically studied in the inner 30 kpc of our galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to 100 kpc. We obtain follow-up spectra of stars in two strong overdensities—including the previously identified outer Virgo Overdensity—and find them to be relatively metal rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming 60–90 kpc counterparts to the 15–20 kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from 50 to 100 kpc, in the same plane as the Sagittarius Stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our galaxy’s stellar halo. 

Abstract Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia–Sausage–Enceladus (GSE), estimated to have occurred ≳8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar overdensities in diagonally opposing octants of the Galaxy (the Hercules–Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several gigayears in the Galactic halo. We employ the simulated stellar halo from Naidu et al., which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many gigayears. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC–VOD. Such halo–disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic center. 

Abstract Several lines of evidence suggest that the Milky Way underwent a major merger at z ∼ 2 with the Gaia-Sausage-Enceladus (GSE) galaxy. Here we use H3 Survey data to argue that GSE entered the Galaxy on a retrograde orbit based on a population of highly retrograde stars with chemistry similar to the largely radial GSE debris. We present the first tailored N -body simulations of the merger. From a grid of ≈500 simulations we find that a GSE with M ⋆ = 5 × 10 8 M ⊙ , M DM = 2 × 10 11 M ⊙ best matches the H3 data. This simulation shows that the retrograde stars are stripped from GSE’s outer disk early in the merger. Despite being selected purely on angular momenta and radial distributions, this simulation reproduces and explains the following phenomena: (i) the triaxial shape of the inner halo, whose major axis is at ≈35° to the plane and connects GSE’s apocenters; (ii) the Hercules-Aquila Cloud and the Virgo Overdensity, which arise due to apocenter pileup; and (iii) the 2 Gyr lag between the quenching of GSE and the truncation of the age distribution of the in situ halo, which tracks the lag between the first and final GSE pericenters. We make the following predictions: (i) the inner halo has a “double-break” density profile with breaks at both ≈15–18 kpc and 30 kpc, coincident with the GSE apocenters; and (ii) the outer halo has retrograde streams awaiting discovery at >30 kpc that contain ≈10% of GSE’s stars. The retrograde (radial) GSE debris originates from its outer (inner) disk—exploiting this trend, we reconstruct the stellar metallicity gradient of GSE (−0.04 ± 0.01 dex r 50 − 1 ). These simulations imply that GSE delivered ≈20% of the Milky Way’s present-day dark matter and ≈50% of its stellar halo. 

Abstract The astrophysical origins of r -process elements remain elusive. Neutron star mergers (NSMs) and special classes of core-collapse supernovae (rCCSNe) are leading candidates. Due to these channels’ distinct characteristic timescales (rCCSNe: prompt, NSMs: delayed), measuring r -process enrichment in galaxies of similar mass but differing star formation durations might prove informative. Two recently discovered disrupted dwarfs in the Milky Way’s stellar halo, Kraken and Gaia-Sausage Enceladus (GSE), afford precisely this opportunity: Both have M ⋆ ≈ 10 8 M ⊙ but differing star formation durations of ≈2 Gyr and ≈3.6 Gyr. Here we present R ≈ 50,000 Magellan/MIKE spectroscopy for 31 stars from these systems, detecting the r -process element Eu in all stars. Stars from both systems have similar [Mg/H] ≈ −1, but Kraken has a median [Eu/Mg] ≈ −0.1 while GSE has an elevated [Eu/Mg] ≈ 0.2. With simple models, we argue NSM enrichment must be delayed by 500–1000 Myr to produce this difference. rCCSNe must also contribute, especially at early epochs, otherwise stars formed during the delay period would be Eu free. In this picture, rCCSNe account for ≈50% of the Eu in Kraken, ≈25% in GSE, and ≈15% in dwarfs with extended star formation durations like Sagittarius. The inferred delay time for NSM enrichment is 10×–100× longer than merger delay times from stellar population synthesis—this is not necessarily surprising because the enrichment delay includes time taken for NSM ejecta to be incorporated into subsequent generations of stars. For example, this may be due to natal kicks that result in r -enriched material deposited far from star-forming gas, which then takes ≈10 8 –10 9 yr to cool in these galaxies.